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ABSTRACT 

A521 is an interacting galaxy cluster located at z=0.247, hosting a low frequency radio halo connected to an 
eastern radio relic. Previous Chandra observations hinted at the presence of an X-ray brightness edge at the 
position of the relic, which may be a shock front. We analyze a deep observation of A521 recently performed 
with XMM-Newton in order to probe the cluster structure up to the outermost regions covered by the radio 
emission. The cluster atmosphere exhibits various brightness and temperature anisotropics. In particular, two 
cluster cores appear to be separated by two cold fronts. We find two shock fronts, one that was suggested by 
Chandra and that is propagating to the east, and another to the southwestern cluster outskirt. The two main 
interacting clusters appear to be separated by a shock heated region, which exhibits a spatial correlation with 
the radio halo. The outer edge of the radio relic coincides spatially with a shock front, suggesting this shock is 
responsible for the generation of cosmic ray electrons in the relic. The propagation direction and Mach number 
of the shock front derived from the gas density jump, M = 2.4 + 0.2, are consistent with expectations from the 
radio spectral index, under the assumption of Fermi I acceleration mechanism. 

Keywords: gs\w6.ts: clusters: general, — galaxies: clusters: individual(A521), — galaxies: clusters: intraclus- 
ter medium, — shock waves 



1. INTRODUCTION 

Collisions between massive galaxy clusters are the most 
energetic events in the present universe. Part of the ki- 
netic energy released during these collisions is dissipated 
through supersonic shock fronts propagating in the intra- 
cluster medium (ICM) and turbulent motions. While heat- 
ing the thermal component of the ICM, shocks and turbu- 
lence may also accelerate (or reaccelerate) relativistic parti- 
cles (e.g. ICassano & Bmnetti 2005; Hoeft & Brugge nll2007t 



Rvu et al.l I2003i: iPfrommer et al.i i2006i: .Bmnetti & Lazariar,, 
20071; IVazza et al.l 120091) . Radio observations probe these 



complex mechanisms through the detection of diffuse syn- 
chrotron emission from the ICM, in the form of giant radio 
halos, Mpc-scale radio emission in the cluster central regions, 
and radio rel ics, sharp-edged ra dio sources in the cluster pe- 
riphery (e.g.JFerrari et al.EOOSi lCassanoll2009t I\^ril20 1 lb 



lBrunettiil201 1 , for recent reviews) 

Observable as sharp X-ray brightness and temperature dis- 
continuities, few shock fronts have been detected so far, be- 
cause they can only be visible in the bright est cluster regions 
and i n parti c ularly favorab l e prqi e ctions (iMarkevitch et alJ 



' "2^05^, 'RusseUeLaU l2()ToE iFinoguenov et alJ 120101: 
2011). Peripheral radio reUcs are believed to 
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be shock fronts that propagated far outside the X-ray bright 
region, while still accelerating (or re-accelerating) electrons, 
which produce radio emission and quickly cool after the 
shock passes, resulting in a characteristic narrow feature (e.g., 
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lEnssUn et allll998l; Ivan Weeren et al.ll2011[) . The physics of 
giant radio halos is probably more complex. Radio halos 
plausibly result from the (re)-acceleration and transport of 
relativistic particles in large turbulent regions of the ICM, 
although many aspects of the mechanisms generating radio- 
emitting electrons remain unclear (e.g., Bru nettill201 1 1 for a 
recent review). Sharp radio edges (and radio relics) are fre- 
quently observed at the border of giant radio halos, suggesting 
a possible link between merger shocks and the generat ion of 
turbulence in the ICM (e.g., ,Markevitch 2010; Macar io et al.l 
l20Tlh . 

A521 is a moderately dis tant (z=0.247) and X-ray luminous 
(Lx = (5.2 -H 1.2) X 10^^W,|^MudeLaL! 2000) □ galaxy clus- 
ter, presenting several signatures of dynamical activity. As 
revealed akeady in ROSAT images, its X-ray and optical com- 
ponents appear spatiall y segregated, with a n N-S bimodality 
of the X-ray emission (lAmaud etal]l2000l) . and a more com- 
plex galaxy number density dist ribution revealing tw o NW/SE 
and NE/SW major elongations (iFerrari et al.ll2003 ). As fur- 
ther shown from Chandra data analysis, the ICM in A521 ex- 
hibits an irregular thermal structure with indications for gas 
heating at the interf ace between the two main gas components 
(IFerrari et al.ll2006l) . A521 exhibits a giant radio halo that is 
the protot ype of the cla ss of ultra-steep spectrum radio ha- 
los (Brunet ti et al.ll20()8l) . The halo being spatially connected 
to a radio relic, A521 provides us with an ideal test case to 
investigate the effects of shocks on the properties of thermal 
and nom-thermal components of the ICM and their connec- 
tion with giant radio halos. In this merging cluster, a shock 
has been suggested by the presence of an X-ray brightness 
edge on the S E side of the cluster, coin ciding with the edge of 
the radio relic (iGiacintucci et al.ll2008l) . The larger-scale clus- 
ter radio halo shows a very steep synchrotron spectrum sup- 
porting a picture where relativistic electrons are stochastically 
re-accelerated by the non-linear interaction with turbulence in 

' 10"W = lO^^erg.s"'; X-ray luminosity in the 0.1-2.4 keV band has 
been corrected for luminosity distance assuming Ho = 70 km s"' Mpc"', 
and A = 0.7. 
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Figure 1. EPIC XMM-Newton exposure of A521. 

the ICM jBrunetti et al.ll2008b . 

The present article will focus on the analysis of a deep ob- 
servation of A521 recently performed with XMM-Newton , 
with particular goal to probe the ICM structure up to the out- 
ermost regions covered by the cluster radio halo, and radio 
relic. After discussing data preparation and analyzes issues 
in Sections |2] and |3] we present the various X-ray brightness 
and temperature features revealed by this observation in Sec- 
tion |4] We comment on the interplay between thermal and 
non-thermal components of the ICM in Section |5] Unless 
otherwise noted, any energy distribution is normalized as a 
probability density function, while confidence ranges on indi- 
vidual parameter estimates are 68 %. In the following, intra- 
cluster distances are computed as angular diameter distances, 
assuming a A-CDM cosmology with Hq = 70 km s"' Mpc"\ 
Qm = 0.3, Qa = 0.7. Given these assumptions, an angular 
separation of 1 arcmin corresponds to a projected intra-cluster 
distance of 232.5 kpc. 

2. OBSERVATIONS AND DATA PREPARATION 

The EPlC-XMM-Newton data set is a dual observation of 
A521, performed with focal aim point and central EPIC-MOS 
CCDs located in the SE cluster outskirts (see Figure[T). In or- 
der to remove the contribution of soft proton flares, we filtered 
the histogram of the photon arrival times through a temporal 
wavelet analysis. A summary of the 'good' exposure time 
remaining on each of the three EPIC cameras is provided in 
Table[T] The average good exposure time is about 75 ks. 

In order to perform imaging and spatially resolved spec- 
troscopy, we binned photons in sky coordinates (k, I) and en- 
ergy (e), matching the angular and spectral resolution of each 
focal instrument. To map the surface brightness of extended 
sources, these photon counts may have to be normalized for 
spatial and spectral variations of the telescope effective area 
and detector exposure times. We thus associate an 'effective 
exposure' array, E{k, I, e), to the photon event cube. Express- 
ible e.g. in s.deg^, E(A;, I, e) is computed as a linear combina- 
tion of CCD exposure times, tccoiK U p), related to individual 
observations p, with local corrections for useful CCD areas. 
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Figure 2. Background spectrum observable in the A521 outskirts. Light 
blue: particle background. Cyan blue: Cosmic X-ray Background emission. 
Blue and Violet: T AE emission (kTl = 0.099 keV, kT2 = 0.248 keV, see 
IKuntz & Snowder] {20001) . and details in Sect. 13.11 1. Green: Residual soft 
proton emission. Red and black: overall fit and data set. 

accoik, l-J>), Reflexion Grating Spectrometer (RGS) trans- 
missionfl frRcs, and mirror vignetting factors accD(^, U /^fl 



Eik, e) = ^ tcCDik, I, p) X Aamirror(^, U e, p) 

X tr^Gsik, l,e,p) 
XaccD{k,l,p) (1) 

3. DATA ANALYSIS 

3.1. Background noise modeling 

The cluster emissivity must be separated from an additive, 
spatially extended and mostly stationary background noise 
including false photon detections due to charged particle- 
induced and out-of-time events, but also the Cosmic X-ray 
Background (CXB), and some Galactic foreground compo- 
nents. 

^ EPIC-MOS detectors share a common optical path with the RGS 
' Information about these instrumental effects have been obtained from the 
XMM-Newton -EPIC Cun-ent Calibration Files (CCFs) 



Shock heating of A521 



3 



Table 1 

Effective exposure time of each XMM-Newton -EPIC observation. 



XMM-Newton Centre coordinates MOSl effective MOS 2 effective PN effective 

obs. IDs exposure time (ks) exposure time (ks) exposure time (ks) 

0603890101 (S) 04h54m22.00s-10"16'30." 15.7 (50.8 %) 15.7 (68.9 %) 12.8 (36.1 %) 

0603890101 (U) 04h54m22.Q0s-10°16'3Q." 64.4(77.3%) 64.5 (74.6 %) 57.6(52.5 %) 



Note. The fraction of the useful exposure time after solar-flare filtering is shown in brackets. 



The XMM EPIC background is dominated by the particle 
component, which is modeled from observations performed 
in the Filter Wheel Closed (FWC) position during revolutions 
230 to 2027 as for the EPIC-MOS cameras, and 355 to 1905 
as for the EPIC-PN camera. Following an approach proposed 
in e.g. Kuntz & Snowden (2008) or Leccardi & MolendJ 
(12008) this model sums a quiescent continuum to a set of flo- 
rescence emission lines convolved with the energy response 
of each detector It is completed with a residual emission as- 
sociated with soft protons, presently only detectable in the 
case of the EPIC-PN camera and modelled as a power-low 
spectrum normalized to 1.4 cts deg"^ s"' in the 0.5-1. keV 
band. To account for two different spectral shapes in the soft 
and hard bands, the quiescent continuum is modelled as the 
product of a power law with an inverted error function in- 
creasing in the soft band. We set the emission line e nergies 
to the values reported in iLeccardi & Molen"d ? ('2008'), while 
the soft proton residual is modelled using an additional power 
law. Presumably due to differences in the collecting areas of 
the imaging and readout detector regions, the EPIC-MOS qui- 
escent continuum exhibit a small emissivity gradient along the 
RAWY CCD coordinate, which has been measured and taken 
into account in the model. Because the fluorescence lines 
exhibit a more complex spatial variation (iLumb et alj 120021 
iKuntz & Snowdenll2008[) . we modelled the emissivity distri- 
bution of the most prominent lines H from the wavelet filter- 
ing of a set of FWC event images in narrow energy intervals 
around each line. 

Secondary background components include the Cosmic X- 
ray background and Galactic foregrounds. Being associated 
with real photon detections, these components are corrected 
for the efifective exposure. The Cosmic X-ray background 
is modell ed with an absor bed power law of index y = 1 .42 
(see, e.g.. Lum b et al.ll200 2). while the Galactic foregrounds 
are modelled by the sum of two absorbed thermal components 
accounting for the Galactic transabsorpti on emission (TAB; 
kT]_ = 0.099 keV and kT2 = 0.248 keV, see lKuntz & Snowded 
120001) . We estimate emissivities of each of these components 
from a "joint-fit" of all background components in a region of 
the field of view located beyond the boundary of X-ray emis- 
sion in the SE cluster outskirt, but covered by the central MOS 
CCDs (see also Figure|2]i. This estimates yields 13.4, 28.0 and 
29.5 cts.m^^deg^^.s"' in the 0.5-1. keV band as for the two 
transabsorption and CXB components, respectively {x^jd.o.f. 
- 1.23). Our background model, ^h<±{k,l)Yjg{k,l,e), even- 
tually includes a contribution for the EPIC-PN out-of-time 
count rate, which is estimated in each energy band as 6.3 % 
of all photon counts registered along the CCD columns. 

3.2. Spectroscopic and surface brightness measurements 

Namely the Al, Si and Cu, Ni complexes as for the EPIC-MOS and 
EPIC-PN cameras, respectively. 



To estimate average ICM temperatures, kT, and metal abun- 
dances, Z, along the line of sight and for a given location of 
the field of view {k, I), we add a source emission spectrum to 
the background model, and fit the spectral shape of the result- 
ing function, Nevt(^, Fevt(kT, Z, Nh, e), to the photon energy 
distribution registered in the energy band (0.3-12 keV): 

Nevt(A:,OFevt(kT,Z,NH,e) 
= E(k, I, e) X Nicm(^, /) FicM(kT, Z, Nh, e) 
+Nbck(A:,/,e). (2) 

In this modeling, the source emission spectrum 
FicM(kT, Z, NH,e) assumes a redshifted and Nh absorbed 
emission modelled from the Astrophysical Plasma Emission 
Code (APEC, Smith et al. 2001), with the element abun- 
dances o flGrevesse & Sau val ( 1998) and neutral hydrogen ab- 
sorption cross sections of Balucinska-Church & McCammo^ 
( Il992l) . The Nh value has been fixed to 4.9 x lO^^m ^ 
from measurements obtained near A521 in t he Lei- 
den/A rgentine/Bonn Survey of Galactic HI (iKalberla et alJ 
120051) . It is altered by the mirror effective areas, filter trans- 
missions and detector quantum efficiency*^, and convolved 
by a local energy response matrix M(k,l,e,e') computed 
from response matrixes files (RMF) tabulated in detector 
coordinates in the XMM-Newton -EPIC calibration data base. 

B 

To compute images and radial profiles of the intra-cluster 
gas distribution, we estimate a cluster surface brightness map, 
2x(^, 0^ from photon counts registered in a given energy band, 
and corrected for effective exposure and additive background. 
Assuming an average ICM energy distribution, < FicM(e) >, 
we define 2x(A:, 1) as a function of an effective exposure map, 
e(^,/) = i:, <FicM(e)>E(^,/,e): 

Yje nevt(^, U e) - Nbck(^, Tje U 
^x(k, I) = . (3) 

e(k,r) 

All parameters of < FicM(e) > = FicM(kTo,Zo,NH,o,e) 
are practically determined from spectral fitting of the main 
cluster emission spectrum: kTo = 6.7 keV, Zo = .4Z0, Nh,o - 
4.9 X 10^'*m"^, while IL^ikJ) is estimated in a "soft" energy 
band ([.5-2.5] keV') in order to lower the dependence of < 
FicM(e) > on kTo. 

3.3. Surface brightness, ICM density and temperature 
profiles 

In the following, surface brightness and temperature pro- 
files have been extracted within cluster sectors oriented ap- 
proximately along the surface brightness gradients. We de- 

' ' EPIC response matrixes are computed from caimed RMFs con'espond- 
ing to the observation period provided by the XMM-Newton Science Opera- 
tion Centre. 
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rived the radial surface brightness £x('") by averaging the sur- 
face brightness Sx(^, of Equation Q in each profile annulus 
composed of N pixels (k, I), as follows: 

IjYj^AkJ), (4) 

kj 

The background contribution Nbck(fc, (k, I, e) being es- 
timated within a much larger area of the field of view than 
any sector annulus used to derive Ilx('"), we neglected any sys- 
tematic uncertainty related to its modeling and estimated the 
variance on 2x(r) from a weighted mean of the local Poisson 
fluctuations in E^ik, 1): 

2 _ 1 y 0'n,Jk,lf _ 1 yi nevt(fc,Q 

Projected temperatures kT(r) and associated confidence in- 
terval (5kT(r) have been computed within each annulus by fit- 
ting a uniform emission model to the data set. To do so, we 
averaged the emission models of Equation associated with 
each pixel (k, I) of the annulus, and estimated the model pa- 
rameters kT(r), Z(r), NH(r), via ax^ minimization. 

These brightness and temperature profiles have been used 
to model the underlying density and temperature of the ICM, 
assuming spherical symmetry of the cluster atmosphere in the 
vicinity of the features of interest. This was undertaken by 
projecting and fitting parametric distributions of the three- 
dimensional (3D) emission measure npWf, and tempera- 
ture, T(r), to the observed profiles. In this modeling, projected 
brightness profiles are convolved with the XMM-Newton fo- 
cus Point Spread Function (PSF), while projected tempera- 
tures are computed assumi ng the 'spectroscop ic-like' weight- 
ing scheme proposed in Mazzo tta et al.l (12004). In Section l431 
the ICM emission measure [3 and temperature profiles across 
two shock fronts have been modelled by step-like distribu- 
tions with a common jump radius rf 

r ^ {X:>lnl {r I r j)-^'^' , r < r j ... 

= | „V(,/,.f2,.^,>,/, (6) 

Tw^rT^;;;^- 

3.4. Imaging and spectral-imaging 
3.4.1. Imaging 

An image of the cluster is presented on the top panels of 
Figure [3] To obtain this image, l.^{k, I), we corrected the 
EPlC-XMM-Newton raw photon image for spatially variable 
effective area and background flux, following Equation (|3). 
The point-like sources have also been modelled by means of 
an isotro pic undecimated B3-spline wavelet analysis (see e.g. 
IStarck et al.. 200 7). and subtracted from the image. 

A map of anisotropic details in the ICM structure is shown 
on the bottom-left panel of Figure[3] To create this image, we 
subtracted a wavelet filtered map of the photon rate, 'Lxik, 0> 
from the photon rate itself, then smoothed the residual image 
with a gaussian function of typical width FWHM = 20 arcsec. 
The wavelet filtering has been performed by means of a soft 
3cr thresholding of B3-spline wavelet coefficients, the sig- 
nificance thresholds being directly computed from the raw 

More precisely, the ICM emission measure per volume unit. 



— Poisson distributed — photon map, followin g the multi- 
scale v ariance stabiUzation scheme introduced in lZhang et alj 
(l2008h . 

3.4.2. Spectral-imaging 

In order to map the ICM temperature in A521, we used 
the EPIC-XMM-Newton data se t and applied the spectral- 
imaging algorithm detaile d in iBourdin et al.l (12004) and 
iBourdin & Mazzottal (12008 1 , hereafter BOS). Following this 
algorithm, a set of temperature arrays kT(A:, /, a) with asso- 
ciated fluctuations <Tit,(k, I, a) are first computed on various 
analysis scales a, then convolved by complementary high- 
pass and low-pass analysis filters in order to derive wavelet 
coefficients. The wavelet coefficients are subsequently thresh- 
olded according to a given confidence level in order to restore 
a de-noised temperature map. Here, the signal analysis have 
been performed over 6 dyadic scales within an angular resolu- 
tion range of 6a - [1.7 - 1 10] arcsec. This was undertaken by 
averaging the emission modelled by Equation (|2]l within over- 
lapping meta-pixels {k, I, a), and computing the kT(A:, /, a) and 
criii(k, I, a) arrays by means of a likelihood maximization. The 
resulting ICM temperature map shown in Figure |3] was then 
obtained from a B2-spline wavelet analysis (see BOS for de- 
tails) with coefficients thresholded to the Icr confidence level. 

4. ICM THERMODYNAMICS 

4. 1 . Intra-cluster gas brightness and thermal structure 

The X-ray photon image of Figure |3]reveals us the complex 
morphology of the intra-cluster gas in A521. On large scales, 
a northern subcluster with comet shape is apparently falling 
on the main component. The photon image also reveals the 
strongly irregular morphology of the surface brightness, pre- 
senting various edges indicated with dashed lines. Some of 
these brightness jumps have been enhanced in the bottom-left 
image of anisotropic details. They are also noticeable on the 
surface brightness profiles of Figure 2]and|5] At the inter- 
face between the two main interacting cluster components, 
we observe two bow-shaped brightness jumps, CFi and CF2, 
joining each other to form a low brightness cross-shaped fea- 
ture. A third brightness jump with higher curvature radius, S 1 , 
is crossing the southern cluster component from SE to NW, 
while a fourth one, S2, is visible at the South-east cluster out- 
skirts. 

The ICM temperature map of Figure |3] is strongly irregu- 
lar, and presents various noticeable features. The northern 
sub-cluster is clearly cool (kT ^ 4.5keV). The interacting 
region separating this cool core from the main cluster to the 
South appears hotter (kT > 7ke V) and strongly disturbed. The 
cross-shaped brightness depression observable on the photon 
image seems to coincide with a hot cross (kT ^ 9keV), in 
particular, along the brightness jump CFi. The southern part 
of the main cluster is cooler (kT ^ 4keV) than the interacting 
region, in particular, to the South of the brightness jump S 1 . 

Bringing together the brightness and temperature maps of 
Figure [3] we observe that the two brightness jumps CFi and 
CF2 are associated with temperature increases as the bright- 
ness decreases, while the brightness jumps S 1 and S2 are as- 
sociated with a temperature decrement. The two jumps CFi 
and CF2 are thus likely to be cold fronts separating the dens- 
est parts of the two sub-clusters from their interacting region, 
while Si and S2 are probably shock fronts propagating out- 
wards from the colliding clusters. 
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Figure 3. EPIC XMM-Newton observation of A521. Top panels: Photon rate image in the .5-2.5 keV band. Photon counts in this image have been corrected 
for spatially variable effective area, background flux and wavelet detected point-like sources. Bottom-left panel: Anisotropic details in the IC M emissivity 
map. These details have been enhanced from subtraction of a wavelet denoised map to the photon rate (further details are provided in Sect. 13.4.11 Bottom-right 
panel: ICM temperature map obtained from wavelet spectral-imaging. Prominent brightness jumps are indicated by dashed Hnes on the photon rate image and 
temperature map. 

One of the most striking features seen in our data is a cross- 
shaped brightness depletion separating the two colliding sub- 
clusters. This feature also corresponds to a temperature and 
entropy enhancement, in particular just outside the two cold 
fronts, CFi and CF2. What we see is probably shocked gas 
with high entropy being squeezed by the converging cool core 
remnants and flowing around the densest part of the two inter- 
acting clusters, without penetrating the two cold fronts. The 
projected layer of shocked gas would thus exhibit maximal 
temperature and entropy near the two cold fronts, where it is 
tangentially intercepted by the line of sight. This shocked gas 
layer might also partly overlay in projection the main sub- 
cluster from its boundary delineated by cold front CF2, to the 
southern brightness jump, S 1 . A possible interpretation for the 
origin of this hot gas flow is illustrated on Figure|6] Originally 
located at the cluster boundary (if there is one), the high en- 



Figure 4. Projected gas brightness measured across two cluster sectors inter- 
cepting the brightness jumps CFi and CFt. The projection of a step-like gas 
density distribution (Equation |6) convolved with the XMM-Newton PSF is 
superimposed as a dotted line, assuming density jump amplitudes of 1 .7 ± . 1 
and 2.1 ± .1 as for CFi and CF2, respectively. 



4.2. 



Cold fronts and shock heating at the interface between 
two interacting sub-clusters 
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Figure 5. Top panels: Projected gas brightness and temperature profiles measured across two cluster sectors intercepting the brightness jumps Si and S2, as 
shown on FigureJT] Bottom panels: ICM density and temperature profiles modelled as step-like 3D distiibutions matching the projected profiles (see also 
Equations [6|and[2f Dispersions on these profiles have been estimated from random realizations of the data set and corresponding models, each profile envelope 
delimiting 68 % of the realizations with closest distance from the original data set. The projection of these distributions is reported as a dotted line on the 
projected profiles. 



Table 2 

Estimated cluster photon counts within the ,3-5, keV band, in the regions shown on Figure|2] 



Detector 


Region 1 


Sector 1 
Region 2 Region 3 


Region 4 


Region 1 


Sector 2 
Region 2 


Region 3 


EPIC-MOSI 
EP1C-MOS2 
EPIC-PN 


985 (79,4%) 
1016 (80,4%) 
2086 (74,0%) 


518(86,1%) 
552 (87,6%) 
1083 (82,3%) 


592 (88,9%) 
578 (89,4%) 
1 146 (86,4%) 


879 (90,4%) 
868 (91,0%) 
1557 (87,3%) 


140(19,6%) 
87 (14,9%) 
192 (13,6%) 


173 (37,0%) 
144 (36,4%) 
284 (29,2%) 


232 (49,8%) 
232 (49,8%) 
420 (42,9%) 



Note. The fraction of the total counts is shown in brackets. 
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Figure 6. Tentative interpretation of the ICM thermal and entropy structure 
observed in the central region of A521, Left: Early stage of a two-cluster 
merger: the cluster boundaiies start to collide and develop two shock fronts 
propagating within the densest regions of each cluster In the meanwhile, the 
cluster develop two cold fronts while pushing the higher entropy gas away 
from their interacting region. Right: The shock fronts have now propagated 
to the most external regions of the interacting clusters, but could not penetrate 
the two cool cores, A shocked gas region with high entropy remains at the 
interface between the two cold fronts. 
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Figure 7. Photon rate image of A521 extracted in the ,5-2,5 keV band. The 
image has been re-binned to a 6,8 arcsec angular resolution in order to en- 
hance the brightness jumps S 1 and S2 , The two annular sectors show the two 
cluster regions where temperature and brightness profiles of Figure [5] have 
been extracted. 



tropy gas may have been shock heated between the two clus- 
ters starting to interact. It would now expand over the clus- 
ter atmosphere, following shock fronts presently propagating 



outside the cluster cores. One of these shock fronts might be 
observed to the South of the main cluster as the brightness and 
temperature jump, S 1 . 
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Figure 8. Left panel: Anis otropic details in the X-ray emissivity (same as Figure |3] overlaid with emissivity iso-contours in the 240 MHz radio band (Giant 
Metrewave Radio Telescope; IBrunetti et aU2008B . Right panel: ICM temperature map overlaid with the same radio emissivity iso-contours as for the left figure. 

Table 3 

Density, temperature jumps and Mach numbers estimated across the shock fronts Sj and St. 



Density estimates Temperature estimates 

Shock front Jump amplitude Mach number Jump amplitude Mach number 

j+nm 1 lo+o.oi 1 7O+0.6J 1 nf.+i.)A'i — 



1.48 
2.64+ 



1.33 

2.42+! 



1.78_^ 
4.47+^ 



1.76 
3.40+ 



L21. 



4.3. Shock fronts propagation to the cluster outskirts 

The 2D gas brightness and temperature maps of Figure |3] 
suggest the brightness jumps Si and S2 to be shock fronts 
propagating outwards the cluster center Located at various 
distances from the cluster center, these two shock fronts might 
have been developed during two successive cluster collisions. 
In order to analyze ICM thermodynamics across these jumps, 
we extracted the brightness and temperature profiles shown 
on Figure |5] corresponding to the two sectors of Figure An 
estimation of the cluster photon counts in each sector region 
is provided in Table|2] The brightness jumps S 1 and S2 exhibit 
the typical shape of a projected spherical density jump, con- 
volved with the XMM-Newton PSF We model the underlying 
gas density and temperature profiles as two step-like functions 
with common jump location, following equations (|6]l and (|7]i. 
A discussion about the validity of the assumption of the ICM 
spherical symmetry in the vicinity of the shock fronts is pro- 
vided in the Appendix. The gas density and temperature dis- 
tributions corresponding to this model are reported under the 
projected profiles on Figure |5] The 3D density and tempera- 
ture jumps associated with these distributions are reported in 
Table [3] with confidence intervals estimated from the 68 % 
percentiles of a parameter sample matching several random 
realizations of the data set. 

The direction of the temperature jumps is consistent with 
the shock front interpretation. The cold front hypothe- 
sis would instead imply a temperature increase across the 
jumps (Dt < 1), which is excluded by the data. Assum- 
ing two shocks propagating outwards in the main cluster, one 



should be able to estimate the shock Mach numbers from the 
Hugoniot-Rankine density, temperature or pressure jump con- 
ditions across the fronts. Such Mach number values are re- 
ported on Table [3] The Mach numbers independently esti- 
mated from the density and temperature jumps are consistent 
with each other, though estimates from the temperature jumps 
have larger uncertainties. We will hereafter use Mach number 
estimates for both shocks S 1 and S2, from their density jumps: 
Msi,p = 1.33!0:02 andMs2,p = 2.42 ± 0.19. 



5. NON-THERMAL ICM EMISSION 

A521 hosts a radio relic in its southeastern peripheral re- 
gion, and a rare low-frequency giant radio halo. In order to in- 
vestigate the interplay between thermal and non-thermal com- 
ponents of the ICM emission, the 240 MHz radio image ob- 
tained from observations performed at the Giant Metrewave 
Radio Telescope (GMRT) has been superimposed on the X- 
ray photon and ICM temperature maps of Figure |8] 

5.1. The A521 radio halo 

The A521 radio halo has been discovered from low- 
frequency observations at the GMRT (240, 325, 610 MHz, 
Brunetti et al. 2008, see also Figure |8]l and then studied in de- 
tail through a deep follow-up Very Large Array observation 
at L4 GHz (Dallacasa et al. 2009). Its very steep spectrum, 
with spectral index a ~ 1.9 between 325 and 1400 MHz, 
suggests magneto-hydrodynamic turbulence to be responsi- 
ble for the in-situ re -acceleration of the relativistic electrons 
dBrunetti et al.l l2008V The radio halo is covering the cluster 
central region, exhibiting an E-W elongation and reaching the 
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Figure 9. Projected galaxy density distribution derived from pliotometric ob- 
servations performed at the CFH telescope (dark color indicate higher den- 
sities, Dressier algorithm, see Ferrari et al. 2003). Pink iso-contours: ra- 
dio emissivity in the 610 Mhz band i Giacintucci et al...2008) . Black iso- 
contours: X-ray emissivity in the .5-2.5 keV band (Curvelet de-noising of the 
XMM-Newton image). 

radio relic to the Southeast. When excluding the relic region, 
the halo appears spatially correlated with the cluster X-ray 
emission. There is an even better correlation between the ra- 
dio brightness and the hottest regions of the ICM, -in partic- 
ular, the radio brightness exhibits a quick drop across the S i 
shock. 

The complex thermodynamics of the ICM in the cluster 
center hint at the possible origin of the turbulence that may 
re-accelerate non-thermal particles in the halo. The two cold 
fronts CFi and CF2 may have developed K-H instabilities at 
large angles from the main cluster collision axis. As sug- 
gested by the spatial correlation between shock-heated re- 
gions and the radio emission, turbulence may alternatively 
have been generated behind the two shocks Si and S2, now 
propagating to the cluster outskirts. In addition, the merger 
disturbance has likely generated turbulence within the two 
subcluster core remnants. 

5.2. Shock front propagation and the radio relic 

A521 has been kno wn to host a SE radio relic observed 
at various f requencies (iFerrari et al.l l2006t IGiacintucci et al.l 
120061120081 hereafter, G08). As shown in G08, the integrated 
synchrotron radiation in the relic exhibit a power-law spec- 
trum with spectral index or ^ 1.5 in the frequency range 235- 
4890 Mhz, with evidences of steepening of the radio spec- 
trum with increasing distance from the eastern edge. As 
further noted in G08, the outer edge of the radio relic co- 
incides with the X-ray edge S2, which we have shown in 
this work to be a shock front propagating to the cluster out- 
skirts. As observed in several peripherical radio sources of 
galaxy clusters (see e.g. Briiggen et al. (2012) for a recent re- 
view), these facts support the shock electron (re)-acceleration 
to be at least partly responsible for the radio emission from 
the relic. Assuming diffuse shock acceleration for the ori- 
gin of the emitting electrons, in the test particle approach the 



slope of the injection spectrum of cosmic rays is related to 
the shock Mach number, M, by dBlandford & Eichleri[l987l) 
6i - 2{M^ + l)/(M^ - 1). This leads to a spectrum of elec- 
trons in the downstream region with slope S - 6i + l (implying 
a synchrotron spectral index a - Sill) taking into account ra- 
diative losses and assuming stationary conditions. G08 thus 
predicted a shock propagation with Mach number, M ^ 2.3, 
from their measurement of a. The steepening of the radio 
spectrum with increasing distance from the eastern edge fur- 
ther allowed them to predict a shock propagation to the cluster 
outskirts. The propagation direction and Mach number of the 
shock front S2, Mp = 2.4 + 0.2 (cf. section |43] ), are fully 
consistent with this hypothesis. 

It is further worth noticing that the X-ray edge correspond- 
ing to the shock front seems to extend in North-South direc- 
tion more than the radio relic. A first interpretation for this 
limited extent of the relic might be that the shock would re- 
accelerate pre-existing relativistic electrons in the ICM. In this 
case the radio relic could reflect the spatial and energy dis- 
tribution of the pre-existing electrons across t he shock front. 
In line with this hypothes is, recent analyzes dKang & Jones! 
120071: iKang & Rvul l2011a 'b') suggest that the presence of pre- 
existing particles in addition to the thermal pool can signifi- 
cantly increase the average efficiency of the particle acceler- 
ation and the expected synchrotron emission at weak shocks 
(M < 3). Differences in extension between the shock and the 
radio relic might alternatively indicate some changes in the 
efficiency of electron acceleration changes along the shock 
front, possibly due to local variations of the Mach number 
(e.g., Hoeft et al.. 2008.) . In this respect the radio relic in 
A521 is located at the extremity of the NW /SE major galaxy 
alignment evidenced in lFerrari et al.l (120031 see also Figure|9]l, 
where indeed recent accretion of subcluster material may have 
produced inhomogeneities in the ICM. 

6. DISCUSSION AND CONCLUSIONS 

A521 is a complex cluster system where optical analyzes 
have revealed at least three galaxy groups to the SE, and four 
groups to the NW including the cluster BCG group coincid- 
ing with the X-ray peak (Ferrari et al. 2003). The X-ray mor- 
phology of the BCG group suggests an infall along a NNW- 
SSE direction (projected onto the sky plane), which is slightly 
offset with respect s to the major NW-SE galaxy alignment 
(IFerrari et al.ll200 6l). The cluster atmosphere exhibits various 
brightness and temperature edges associated with cold fronts 
and shock fronts, that our XMM-Newton data revealed. 

The main two interacting gas components in the central 
region of this system are separated by a region of gas with 
lower density, higher temperature and entropy. We interpret 
this feature a flow of high-entropy gas being squeezed by 
two converging subcluster cores that are delimited by cold 
fronts. We suggest this high entropy gas to have been heated 
by shocks formerly developed when the two gas components 
started to interact. One of these shocks is currently observed 
to the South of the main component, with Mach number 
Mp = 1.33^QQ^. The hot gas region separating the two inter- 
acting components appears spatially correlated with the clus- 
ter radio halo. The development of turbulence in the hot gas 
flowing between the two cool cores may be responsible for 
high energy electron re-acceleration, yielding the radio halo 
emission. Merger shock propagation and/or cold fronts may 
have contributed to the development of these instabilities. 

A shock front is observed at the Southeast cluster outskirt. 
An X-ray brightness edge there has been hinted at by Chandra 
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data (G08), though the statistical significance was marginal. 
The orientation of this shock front and its large distance from 
the cluster center suggest that it is associated with a cluster 
collision that has occurred prior to the current two-component 
interaction. Our Mach number for this shock, Mp = 2.4 + 0.2, 
is consistent with that expected from the spectrum of the ra- 
dio relic in G08 under the assumption of Fermi I accelera- 
tion mechanism. As observe d in X-ray follo w-ups of other 
radio relics -A3667, Finoguen ov et all fcolOb : RXCJ1314.4- 
2515. Mazzotta et aP (l201 1)-. its detection supports the shock 
electron (re)-acceleration to be at least partly responsible for 
the radio emission from the relic. The detection of a polar- 
ization of the relic would be an additional support for this 
process, complementary to the extension of its synchrotron 
spectrum to very high radio frequencies, and to evidences for 
spectral steepening downstream to the shock (G08). Delim- 
ited by the shock front, the radio relic seems however to sub- 
tend only a fraction of the shock front. Differences in the 
spatial extent of a radio relic and its companion shock front 
have a lso been observed in t he colliding cluster RXCJ 1 3 14.4- 
2515 ( IMazzotta et al.ll20Tl1) . where a radio relic seems to be 
confined to a small section of the shock front presumably 
distorted by a nonuniform gas flow. These differences may 
thus reflect variations of the efficiency of particle acceleration 
across the shock that could be driven by local variations of the 
Mach number and shock velocity. Deeper X-ray or SZ obser- 
vations may enable us to investigate this hypothesis, though 
the present XMM-Newton image does not seem to evidence 
any strong variation in the amplitude of the surface bright- 
ness edge, and thus in the shock Mach number An alternative 
hypothesis is that the radio relic would reveal us local inho- 
mogeneities in the properties of the pre-existing relativistic 
electrons, that would be re-accelerated by the shock passage. 
The observed connection between the radio halo and the relic 
may suggests that pre-existing relativistic electrons have first 
been accelerated by turbulent gas motions responsible for the 
radio halo emission, then re-accelerated at the shock front. 
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APPENDIX 

The ICM density and temperature distributions i nterc epting 
the shocks Si and S2 has been modelled in Section [33] as two 
step-like functions, assuming the shock center of curvature 
and the ICM centroid coincide to coincide with each other 

The X-ray image of Figure |3] seems however to show that 
the shock fronts Si and S2 are less curved than the closest 



cluster brightness isophotes. To investigate the systematic un- 
certainties inherent to our spherical symmetry approximation, 
we alternatively tried to model the shock front and the ICM 
density as two spherical distributions with distinct centers. 
Assuming these two centers to be located in the plane of the 
sky, the ICM emission measure is now expressed per volume 
unit, as: 



[npne](r) = 




where r and r' refers to the norm of each radius vector in 
the shock and ICM frames, respectively. Introducing d and 
d', the projection of these radius vectors onto the sky plane, 
a surface brightness profile intercepting the shock is obtained 
from integration of Equation ([T]i along the line of sight: 

S(x) = 2x J J J A(T(d,l))[npne](d,d',l,(l>)dddld(l), 

<t>,„i„ x-Sx 

(2) 

where d' is related to c/ as a function of r„, the distance sep- 
arating the shock from the center of the ICM distribution and 
0, the angle separating the projected radius vector to the shock 
propagation axis {d' - ^J(rj - r„)' - 2(rj - ro)d cos(0) + d^, 
see also Fig. [TOb . In addition to an ICM density slope, v, and 
the shock curvature radius, density and temperature jumps, 
rj, D„ and Dt, respectively, the ICM emission measure thus 
depends on an asphericity parameter: 1 - rg/rj. 

We tried to invert [n^ne](r) and its parameters from a mini- 
mization of the distance separating (Equation|2]) from 
the X-ray surface brightness profiles extracted across each 
shock front (see Figure |5). Some of the searched parame- 
ters being degenerated with one another, we first performed 
this inversion by fixing the asphericity parameter to and 
0.5, corresponding to shocks located at distances of rj and 
2 X rj from the cluster center, respectively. We subsequently 
left all parameters free to vary and report the results of our 
measurements in Table HI the confidence interval on each pa- 
rameter being estimated from the 68 % percentiles of a pa- 
rameter sample matching several random realizations of the 
data set. As expected, the shock curvature radius, density and 
temperature jumps obtained when fixing the asphericity to 
are consistent with their estimates derived from the spheri- 
cal model of section [33] A marginal difference in the ampli- 
tude of the density jump is still noticeable, since Equation ([T]l 
yields D„ - 1. 39^'^'^* as for shock Si, while Equation ^ 
yields D„ = 1.33^QQg. This difference is probably related 
to the lack of any variation of the ICM density slope at the 
shock crossing, following Equation ([T]i. Fixing the asperic- 
ity to 0.5 instead of also marginally affect the density jump, 
essentially due to the degeneracy between the ICM aspheric- 
ity and density slope. This degeneracy is noticeable in the 
case of S 1 , the shock front observed with the highest statistics. 
Leaving the ICM asphericity free to vary yields estimates of 
0.45 and 0.75 in the case of Si and S2, respectively, consistent 
with the shock curvature radii observed on the X-ray image 
of Figure |3] The shock density, temperature jump and Mach 
numbers derived from these various assumptions are in any 
case consistent with one another, and wit h the ir estimates ob- 
tained from the sperical model of section ITJI Given the lim- 
ited statistics available, it is difficult to break the degeneracy 
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V 



Shock frame (xyz) pai'ameters 
ij Shock curvature radius 
r Radius vector 
d Projected radius vector 
1 Distance to the sky plane 

Projected angle to the propagation axis 



ICM frame (x'y'z') parameters 
To Shock distance to the cluster centre 
r' Radius vector 
d' Projected radius vector 



Figure 10. ICM geometry across the shocks. Left: Cluster volume cut along the line of sight. Right: Cluster volume cut within the sky plane. 



Table 4 

ICM asphericity parameters across the shock fronts SI and S2. 
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between the ICM asphericity, ICM density slope and shock 
curvature radius in the vicinity of the shocks. For simplic- 
ity purposes, we consequently adopted the spherical model of 
section 13.31 in order to derive the amplitudes of the density 
jumps and Mach numbers of the two shocks Si and 82- 
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